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Abstract — We revisit a new type of a Voronoi diagram, in 
which distance is measured from a point to a pair of points. 
We consider a few more such distance functions, based on 
geometric primitives, and analyze the structure and complexity 
of the nearest- and furthest-neighbor Voronoi diagrams of a 
point set with respect to these distance functions. 

Keywords -distance function; lower envelope; Davenport- 
Schinzel theory; crossing-number lemma 

I. Introduction 

The Voronoi diagram is one of the most fundamental 
concepts in computational geometry, which has plenty of 
applications in science and industry. Much information in 
this respect can be found in |4| and [13]; for important recent 
achievements, see (8). 

The basic definition of the Voronoi diagram applies to a 
set S of n points (also called sites) in the plane: its nearest- 
neighbor Voronoi diagram V(S) is a partition of the plane 
into n regions, each corresponding to a distinct site s € S, 
and consisting of all the points being closer to s than to any 
other site from S. Similarly, the furthest-neighbor Voronoi 
diagram of S is obtained by assigning each point in the plane 
to the region of the most remote site. These notions can be 
generalized to higher-dimensional spaces, different types of 
sites, and in other ways. 

One of the recent generalizations of this concept is a 
family of so-called 2-site Voronoi diagrams Q, which are 
based on distance functions that define a distance from a 
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point in the plane to a pair of sites from a given set S. 
Consequently, each Voronoi region corresponds to an (un- 
ordered) pair of sites from S. The original motivation for the 
study |5| was the famous Heilbronn's triangle problem lfl4l . 
Other motivations are mentioned therein. 

For S being a set of points, Voronoi diagrams under a 
number of 2-site distance functions have been investigated, 
which include arithmetic combinations of point-to-point dis- 
tances 0, ifTTl and certain geometric distance functions 0, 
0, 0. In this work, we develop further the latter direction. 

Let S C M 2 , and consider p,q G S and a point v in the 
plane. We shall focus our attention on a few circle-based 
distance functions: 

• radius of circumscribing circle: C(v,(p,q)) = 
Rad(o(ii,p ; q)), where o(v,p,q) is the circle defined 
by v,p, q and Rad(c) is the radius of the circle c; 

• radius of containing circle: K(v,{p,q)) = 
Rad(C(v,p, q)), where C(v,p,q) is the minimum 
circle containing v,p, 

• view angle: V(v, (p,q)) — £pvq, or, equivalently, half 
of the angular measure of the arc of o(v,p,q) that the 
angle /.pvq subtends; 

• radius of inscribed circle: 1Z(v, (p, q)) is the radius of 
the circle inscribed in A(v,p, q); 

< center-of-circumscribing-circle-based functions: let 
o vpq denote the center of the circle o(v,p,q); then 

o o o 

S(v,(p,q)), A(v,(p,q)), and V(v,(p,q)) are the 
distance from o vpq to the segment pq, the area of 
Ao vpq pq, and the perimeter of Ao vpq pq, respectively; 
and on a parameterized perimeter distance function: 

1 Obviously, o(v,p,q) 7^ C(v, p, q) if any of the three points is properly 
contained in the circle whose diameter is defined by the two other points. 



• parameterized perimeter: V c (v, (p,q)) = \vp\ + \vq\ + 
c ■ \pq\, where c > —1. 

The first and third circle-based distance functions were 
first mentioned in |[T0l . The last function generalizes the 
perimeter distance function V(v,(p,q)) = Per(A(w,p, q) 
introduced in J5], and later addressed in Q, ||9]- 

Since two points define a segment, any 2-point site dis- 
tance function d(v, (p,q)) provides a distance between the 
point v and the segment pq, and vice versa. Consequently, 
geometric structures akin to 2-site Voronoi diagrams can 
arise as Voronoi diagrams of segments. This alternative 
approach was independently undertaken by Asano et al., 
and the "view angle" and "radius of circumscribing circle" 
distance functions reappeared in their works [0, flS] on 
Voronoi diagrams for segments soon after they had been 
proposed by Hodorkovsky ifTUl in the context of 2-site 
Voronoi diagrams. However, as Asano's et al. research was 
originally motivated by mesh generation and improvement 
tasks, they were mostly interested in sets of segments repre- 
senting edges of a simple polygon, and thus, non-intersecting 
(except, possibly, at the endpoints), what significantly alters 
the essence of the problem. 

In this paper, we analyze the structure and complexity 
of 2-site Voronoi diagrams under the distance functions 
listed above. Our obtained results are mostly of theoretical 
interest. The method used to derive an upper bound on the 
complexity of the nearest-neighbor 2-site Voronoi diagram 
under the "parameterized perimeter" distance function is first 
developed for the case of c = 1, yielding a much simpler 
proof for the "perimeter" function than the one developed 
in |9), and then generalized to any c > 0. We summarize 
our new results in Table [] 

Throughout the paper we use the notation Vf 1 (S) (resp., 
Vj: (S)) for denoting the nearest- (resp., furthest-) 2-site 
Voronoi diagram, under the distance function T, of a point 
set S. The set S is always assumed to contain n points. 

II. Circumscribing Circle 

Let o{p,q,r) denote the unique circle defined by three 
distinct points p, q, and r in the plane. We now define the 
2-site circumscribing-circle distance function: 

Definition 1: Given two points p, q in the plane, the 
"circumcircle distance" C from a point v in the plane 
to the unordered pair (p,q) is defined as C(v,(p,q)) = 
Rad(o(«,p,g)). 

For a fixed pair of points p and q, the curve C(v, (p, q)) = 00 
is the line pq. This implies that all the points on pq belong 
to the region of (p, q) in {S). In this section we assume 
that the points in S are in general position, i.e., there are 
no three collinear points, and no three pairs of points define 
three distinct lines that intersect at one point. The given sites 
are singular points, that is, for any two sites p, q, the function 
C(v, (p, q)) is not defined at v — p or v = q. 



Theorem 1: Let S be a set of n points in the plane. The 
combinatorial complexity of • (S) is f2(n 4 ). 

Proof: The n points of S define 0(n 2 ) lines, which 
always have 9(n 4 ) intersection points. All these intersection 
points are features of (S), and hence the lower bound. 

■ 

Theorem 2: Let S be a set of n points in the plane. The 
combinatorial complexity of both V^ n \S) and V^\S) is 
0(n i+e ) (for any e > 0). 

Proof: Clearly, the combinatorial complexity of 
V^ n) (S) or v£ f) (S) is identical to that of the respective 
diagram of the 2-site distance function C 2 (v,(p, q)) = 
Rad 2 (o(u,p, q)). It is known that Rad 2 (o(u,p, q)) = 
((\vp\\vq\\pq\)/(4\Avpq\)) 2 = (((v x - p x ) 2 + (v y - 
Py) 2 )((vx ~ q x ) 2 + (v y - q y ) 2 ){{Px - q x ) 2 + {v y ~ 
qy) 2 ))/(4(vx(py - q y ) - p x {v y - q y ) + q x {v y - p v )) 2 ). 
The respective collection of Q(n 2 ) Voronoi surfaces fulfills 
Assumptions 7.1 of 1 16 p. 188]: 

1) Each surface is an algebraic surface of maximum 
constant degree; 

2) Each surface is totally defined (this is stronger than 
needed); and 

3) Each triple of surfaces intersects in at most a constant 
number of points. 

Hence, we may apply Theorem 7.7 of [ibid., p. 191] and 
obtain the claimed bound on the complexity of (S). 

■ 

III. Containing Circle 

Let C(p, q, r) denote the minimum-radius circle contain- 
ing three points p, q, and r in the plane. (That it, C(p, q, r) 
is the minimum circle containing the triangle Apqr.) We 
now define the 2-site containing-circle distance function: 

Definition 2: Given two points p, q in the plane, the 
"containing-circle distance" K. from a point v in the plane 
to the unordered pair (p,q) is defined as JC(v, (p,q)) = 
Rad{0(v,p,q)). 

In our context we have that p ^ q. Assume first that 
v 7^ p, q. Observe that if all angles of Apqr are acute 
(or Apqr is right-angled), then C (p, q, r) is identical to 
o(p, q, r). Otherwise, if one of the angles of Apqr is obtuse, 
then C(p, q, r) is the circle whose diameter is the longest 
edge of Apqr, that is, the edge opposite to the obtuse angle. 
If v coincides with either p or q, then C(v,p, q) is the circle 
whose diameter is the line segment pq. 

Theorem 3: Let S be a set of n points in the plane. The 
combinatorial complexity of Vf?' (S) is Q(n). 

Proof: For simplicity assume that each point from S 
has a unique closest neighbor in S. For each point p <E S, 
consider its closest neighbor q. Then, the points on the line 
segment pq lying sufficiently close to p belong to the region 
of (p, q) in V^"' ) (S'), which is thus non-empty. Since no 
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Table I 

Our results: Worst-case combinatorial complexities of 2-site Voronoi diagrams of a set S of n points with respect to 

DIFFERENT DISTANCE FUNCTIONS 




(c) Obtuse q 



Figure 1. If p, q have a non-empty region in Vv- (S), then pq is an edge 
in DT(S). 

region is thereby encountered more than twice, vj^\S) has 
at least |~n/2] non-empty regions. The claim follows. ■ 
Theorem 4: Let 5 be a set of n points in the plane. The 
combinatorial complexity of (S) is 0{n 2+£ ) (for any 

£ > 0). 

Proof: Let a point v belong to a non-empty region of 
(p, q). No matter if the triangle Avpq is acute (Figure [TJ a)), 
Avpq is obtuse with v being the obtuse vertex (Figure Qlb)), 
or Avpq is obtuse with p or q being the obtuse vertex 
(Figure [He)), the circle C(v,p,q) cannot contain any other 
point x £ S. Otherwise, regardless of the location of x in 
C(v,p,q), we will always have /C(v, (p,q)) > IC(v, (x,q)), 
which is a contradiction. This follows from the fact (see J6] 
Lemma 4.14]) that given a point set K and its minimum 
enclosing circle C, where G is defined by three points 
a,b,c £ K (resp., two diametrical points s, t £ K), 
removing from K one of a,b,c (resp., one of s, t) will 
result in a point set with a smaller minimum enclosing circle. 
Thus, there is a circle containing p, q that is empty of any 
other site from S. This immediately implies that pq is an 
edge of the Delaunay triangulation of S. Consequently, there 
are 0(n) pairs of sites in S that have non-empty regions 



in Vf? 1 (S). Furthermore, it follows from the definition of 
)C(v, (p, q)) that the respective Voronoi surface of (p,q) 
is made of a constant number of patches, each of which 
is a "well-behaved" function in the sense discussed in 
the proof of Theorem |2] Again, by standard Davenport- 
Schinzel machinery, the combinatorial complexity of the 
lower envelope of these 0(n) surfaces is 0(n 2+e ) (for any 
e > 0), and the claim follows. ■ 
Theorem 5: Let S be a set of n points in the plane. The 
combinatorial complexity of V^f'(S) is 0(n 4+£ ) (for any 
e > 0). 

Proof: As in the proof of Theorem [2] we prove this 
claim by using the upper envelope of Q(n 2 ) "well-behaved" 
Voronoi surfaces. ■ 

IV. View Angle 

We now define the 2-site view-angle distance function: 
Definition 3: Given two points p, q in the plane, the 
"view-angle distance" V from a point v in the plane to the 
unordered pair (p, q) is defined as V(v, (p, q)) = Z.pvq. 
Similarly to the circumcircle-radius distance function, the 
view-angle function is undefined at the n given points. For 
a fixed pair of points p and q, the curve V(v, (p, q)) = tt is 
the open line segment connecting the two points p and q, 
while the curve V(v, (p, q)) = is the line pq excluding 
the closed line segment pq. The curve V(v, (p, q)) = ir/2 is 
the circle whose diameter is the line segment pq (excluding, 
again, p and q). 

Theorem 6: Let S be a set of n points in the plane. The 
combinatorial complexity of V^ n \s) is r2(n 4 ). 

Proof: Consider a set S of n points in the plane. 
An example of the intersection of the complements of two 
segments defined by two pairs of points (with respect to the 
supporting lines) is shown in Figure a). These intersection 
points are features of Vy n \s); we show that there are f2(n 4 ) 
such points. To this aim we create a geometric graph G 
whose vertices are the given points, in which each segment's 
complement defines two edges. We add one additional point 
far away from the convex hull of S, and connect it (without 
adding intersections) to all the rays as shown in Figure |2b). 
We can now use the crossing-number lemma for bounding 
from below the number of intersections of the original rays. 
The lemma tells us that every drawing of a graph with 




(a) Intersection point 




(b) Graph 

Figure 2. The graph of V^ n) (S). 



n vertices and m > An edges (without self or parallel 
edges) has il(m 3 /n 2 ) crossing points |T), lfl2l . In our case 
m = 2Q) = n(n — 1), so the number of intersection points 
in G is Cl(n e /n 2 ) = f2(n 4 ). All these intersection points are 
features of Vy (S), and hence the lower bound. ■ 

Theorem 7: Let S be a set of n points in the plane. The 
combinatorial complexity of both Vy(S) and Vy\s) is 
0{n A+£ ) (for any e > 0). 

Proof: For analyzing V^ n) (S) and V-f^iS) we con- 
sider the function ( — cos Lpvq) instead of that of Lpvq. 
This is permissible since the cosine function is strictly 
decreasing in the range [0,7r], By the cosine law, we have 
— cos /.pvq = {\pq\ 2 — \vp\ 2 ~ \vq\ 2 )/ (2\vp\\vq\). As we have 
already seen more than once in this paper, this means that 
the respective collection of 0(n 2 ) Voronoi surfaces fulfills 
Assumptions 7.1 of |[T6l p. 188]. Hence, we may apply 
Theorem 7.7 of [ibid., p. 191] and obtain the claimed bound 
on the complexity of Vy^'(S). ■ 

Theorem 8: Let S be a set of n points in the plane. The 

( f\ 

combinatorial complexity of Vy (S) is Q(n ). 

Proof: Given a set S of n points in the plane, we 
count the intersections of pairs of line segments, where 
each segment is defined by points of S (see Figure HJa)). 
We create a geometric graph whose vertices are the given 
points, and the edges are the line segments connecting 
every pair of points (see Figure |3jb)). The intersections of 
the segments defined by all pairs of points define features 
of Vy (S), because along these segments the view-angle 
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(a) Intersection point 
1 2 




(b) Graph 

Figure 3. The graph of V^'(S). 




Figure 4. TZ(v, (p, q)) is the radius of the circle inscribed in Avpq. 

function assumes its maximum possible value, ir. We can 
now use the crossing-number lemma for counting these 
intersections. The graph with n vertices and m > An edges 
(without self or parallel edges) has Sl(m, 3 /n 2 ) crossing 
points [fl~), |[T2l . In this case m = Q) = n(n — l)/2, hence 
Q(n 4 ) is a lower bound on the complexity of Vy (S). ■ 
Results by Asano et al. J2) immediately imply that the 
edges of Vy^(S) represent pieces of polynomial curves 
of degree at most three. However, the structure of the part 
of \4 f) (S) that lies outside the convex hull CU{S) of S 
is fairly simple: it is given by the arrangement of lines 
supporting the edges of CH(S). This arrangement can be 
computed by a standard incremental algorithm in optimal 
0(fc 2 ) time and space, where k denotes the number of 
vertices of CH(S). Each cell of the arrangement should then 
be labeled with a pair of sites from S, to the Voronoi region 
of which it belongs; this extra task can be completed within 
the same complexity bounds. 

V. Radius of Inscribed Circle 

We now define the 2-site "radius-of-inscribed-circle" dis- 
tance function: 

Definition 4: Given two points p, q in the plane, the 
"inscribed radius distance" 1Z from a point v in the plane to 
the unordered pair (p, q), denoted by lZ(v, (p, q)), is defined 
as the radius of the circle inscribed in the triangle Avpq 
(Figure 0). 



Theorem 9: Let S be a set of n points in the plane. The 
combinatorial complexity of (S) is il(n 4 ). 

Proof: The intersection point of any two lines defined 
by the points from S is a distinct feature of the Voronoi 
diagram under discussion. Thus, n points in S define 9(n 2 ) 
lines, which have 9(n 4 ) intersection points. ■ 

Theorem 10: Let S be a set of n points in the plane. The 
combinatorial complexity of both (S) and vf^ (S) is 
0(n A+£ ) (for any e > 0). 

Proof: Let p, q be two points in S, and u a 
point in the plane. It is a well-known fact that 
K(v,{p,q)) = 2A{v,(p,q))/V(v,(p,q)), where 
A(v,(p,q)) and V(v,(p,q)) are the area and perimeter, 
respectively, of the triangle Avpq. Both the numerator 
and denominator of this fraction can be written as 
algebraic expressions using the coordinates of the points 
v,p,q. Hence, as above, the standard Davenport-Schinzel 
machinery can be applied for obtaining the claim bounds. 

■ 

Theorem 11: Let S be a set of n points in the plane. The 
combinatorial complexity of (S) is Q(n) in the worst 
case. 

Proof: The complexity of vj^ (S) can be as high as 
Q(n). Let S be a set of n point in convex position with no 
three collinear points. Let p and q be two antipodal vertices 
of CH(S), the convex hull of S, and consider two parallel 
lines £ p 9 p and £ q 3 q tangent to CH(S) only at p and 
q, respectively. Next, consider any point v £ £ p , and let it 
move along £ p in either direction. In the limit, the distance 
from v to any pair (s, t) of sites in S equals the width of the 
infinite strip bounded by two lines parallel to £ p and passing 
through s and t, respectively. Consequently, the points of £ p 
lying sufficiently far from p belong to the Voronoi region 
of (p, q). Since the number of pairs of antipodal vertices of 
CU(S) is 9(n), the bound follows. ■ 
A similar reasoning leads to a conclusion that (S) 
has at most a linear number of unbounded regions. To 
demonstrate this, consider any point u in the plane, and a 
line £ 3 u. Observe that the points of £ lying sufficiently far 
from u belong to the Voronoi region of the pair(s) of points 
from S that define the width of S in the direction orthogonal 
to £, and, thus, represent a pair (pairs) of antipodal vertices 
of CH(S). Since the union of all such lines gives the whole 
plane, and the number of antipodal vertices of CH(S) is at 
most linear, the claim follows. 

VI. Distances Based on the Center of the 
Circumscribing Circle 

Let v,p, q be three points in the plane. Consider the circle 
o(v,p,q) passing through v,p,q with center o vpq . We now 
define three more distance functions based on the above 
notation: 

Definition 5: Given two points p, q in the plane, the 

o o 

three distances, denoted by S(v,(p,q)), A(v,(p,q)), and 




Figure 5. The circle o(v, p, q) is defined by the points v, p, q, and has the 

o 

center at o vpq . S(v , (p, q)) is the distance from o vpq to the segment pq (or, 

o 

equivalently, the height of Ao vpq pq perpendicular to pq), and A.{v, (p, q)), 

o 

and "P{v, (p,q)) are the area and the perimeter of Ao vpq pq, respectively. 

o 

V{v, (p, <?)), respectively, are the distance from o vpq to the 
line segment pq, the area of the triangle Ao vpq pq, and the 
perimeter of Ao vpq pq, respectively (Figure |5}. 

The upper bound of 0(n 4+£ ) (for any e > 0) on the 
complexity of the nearest- and furthest-neighbor Voronoi 
diagrams under each of these distance functions can be, 
again, derived by means of Davenport-Schinzel machinery. 
Below we provide some lower bounds. First, we address the 
nearest-neighbor case. 

Theorem 12: Let S be a set of n points in the plane. The 
combinatorial complexity of (S) and vi (S) is f2(n 4 ) 

5 A 

in the worst case. 

Proof: The key observation is the following. Consider 
a pair (p,q) of sites, and let o(p,q) denote the circle with 
the diameter pq. Then, for any point v 6 o{p, q) \ {p, q}, we 

o o 

have S(v, (p, q)) = A{v, (p, q)) = 0. 

Consider two parallel lines and l-z, and let d denote the 
distance between them. For a given n > 2, let us construct 
a set S of n points as a union of two sets Si C h and 
Sa C h consisting of \n/2] and [n/2\ points, respectively, 
in the following way. The sets Si and S2 are constructed 
iteratively; at each odd step, a new point is added to Si, 
and at each even one — to S2. For any i: 2 < i < n, let 
S\ and S2 denote the two sets constructed so far, and let 
M 1 = {o(p,q)\p 6 S\,q £ S^} denote the set of circles 
defined by pairs of points from different sets. We want 
each circle from M n to pass through precisely two points 
from S (those defining it), each two circles from M n to 
intersect, and no three of them to pass through the same 
point not contained in S. Then 0(n 2 ) circles composing 
M n will give rise to Q(n 4 ) distinct intersection points, each 
belonging to a separate feature of either Voronoi diagram 
under consideration, and the claim will follow. 

To ensure the first property, we select the points so that 
the distance between each two points contained in the same 
set Si is much smaller than d, where i = 1, 2. To guarantee 
the second property, at each step j: 3 < j < n, when adding 
a new point s to the respective set, we make sure that for any 
point t from the other set, the circle o(s, t) passes neither 



through any point from S^ -1 U \ {t} nor through any 

intersection point of the circles from M J_1 . This completes 

the proof. ■ 

Theorem 13: Let S be a set of n points in the plane. The 

combinatorial complexity of vf (S) is f2(n) in the worst 

v 

case. 

Proof: A linear lower bound in the worst case for 

Vo (S) can be obtained in the following way. Choose the 

V 

set S of points to lie on some line I, so that the distance 
between any two consecutive points is 1 . Then, the minimum 

o 

possible value for the distance function J> is obviously 2, and 
can be achieved only for a pair (p, q) of consecutive points. 
For each such pair (p, q), consider the circle o(p, q) with the 
diameter pq. Evidently, for any point v G o(p, q) \{p, q}, we 

o 

have V(v, (p, Q)) = 2, and for any other pair (s, t) of sites, 

o 

V(v, (s,t)) > 2. We conclude that each pair of consecutive 

points along £ has a non-empty region in V„ (S). Since 

V 

there are n— 1 pairs of consecutive points, the bound follows. 

■ 

Second, we address the furthest-neighbor Voronoi dia- 
grams. 

Theorem 14: Let S be a set of n points in the plane. 
The combinatorial complexity of all of V ( J ] (S), V ( J ] (S), 

,n S ' A ' 

and V ( J\S) is ft(n 4 ). 
v 

In each case, the proof is identical to that of Theorem |9] 

VII. Parameterized Perimeter 

Finally, we define the 2-site parameterized perimeter 
distance function: 

Definition 6: Given two points p, q in the plane and a real 
constant c > — 1, the "parameterized perimeter distance" V c 
from a point v in the plane to the unordered pair (p, q) is 
defined as V c (v, (p,q)) = \vp\ + \vq\ + c ■ \pq\. 

We require that c be greater than or equal to —1 since 
allowing c < — 1 would result in negative distances. Letting 
c = — 1 results in a distance function that equals for 
all the points on the line segment pq. If c = 0, we deal 
with the "sum of distances" distance function introduced 
in and recently revisited in lfl7l . For c = 1, the 
above definition yields the "perimeter" distance function 
V(v, (p,q)) = Pcr(Avpq). 

In 19) it was proven that the combinatorial complexity 
of the nearest-neighbor 2-site perimeter Voronoi diagram 
of a set of ?i points is slightly superquadratic in n. In a 
nutshell, the proof was based on the observation that any 
pair of sites that has a non-empty region in the perimeter 
diagram also has a non-empty region in the sum-of-distances 
diagram. This immediately implies that the number of such 
pairs is linear in n. (However, unlike in the sum-of-distances 
diagram, a region in the perimeter diagram is not necessarily 
continuous. We were able to construct examples in which 
the number of connected components of a single region is 
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Figure 6. An empty circle containing sites in V. 

comparable to the number of points!) Again, one can apply 
the standard Davenport-Schinzel machinery and conclude 
the claimed upper bound on the complexity of the diagram. 
It remains unclear whether the worst-case complexity of 
the diagram is linear, quadratic, or in between. The proof 
in (9) of the main observation was extremely complex. We 
provide here an alternative and much simpler proof of the 
same bound, which generalizes to the case of "parameterized 
perimeter" distance function for any c > 0. 

Theorem 15: Let S be a set of n points in the plane. The 
combinatorial complexity of Vp(S) is 0(n 2+£ ) (for any 
e > 0). 

Proof: Refer to Figure [6] Let p, q <E S be two sites 
which have a non-empty region in Vp l \S), and let v be a 
point in this region, noncollinear with p and q. In addition, 
let £ be the perpendicular bisector of the line segment pq. 
Assume, without loss of generality, that \vp\ < \vq\. 

Consider the ellipse O vpq passing through q with v and p 
as foci. By definition, for any point s inside this ellipse we 
have I us I + \ps\ < \vq\ + \pq\. Therefore, 

V(v,(p,s)) = \vs\ + \ps\ + \vp\ (1) 
< \vq\ + \pq\ + \vp\ = V(v,(p,q)). 

This means that s cannot be a site in S, for otherwise v 
would belong to the region of (p,s) instead of to the region 
of (p,q). It follows that the ellipse O vpq is empty of any 
sites other than p and q. 

Now consider the line £' that is tangent to O vpq at q, and 
the ray r perpendicular to £' at q and passing through O vpq . 
It is a known property of ellipses that this ray bisects the 
angle £vqp, and, thus, it intersects the line segment vp, 
say, at point o. The circle C centered at o and passing 
through q is tangent to O vpq at q (as well as at another 
point), and is entirely contained in O vpq . Since v is closer 
to p than to q (by our assumption), it follows that the 
circle C also contains p. (If p were on the extension of 
vp in the shaded area, a contradiction would easily be 
obtained by using the triangle inequality: \op\ > \oq\, hence 



\vp\ = \ov\ + \op\ > \ov\ + \oq\ > \vq\, contradicting the 
assumption that \vp\ < \vq\.) Since O vpq is empty of sites 
(except p and q), so is the circle C. Therefore, pq is an 
edge of the Delaunay triangulation of S. The number of 
such edges is linear in n, the cardinality of S. 

Hence, there are O(n) respective surfaces of these pairs of 
sites. One can now apply the standard Davenport-Schinzel 
machinery (as in the proof of Theorem|2). The claim follows. 

■ 

Finally, we state the following theorem. 

Theorem 16: Let S be a set of n points in the plane. 

(a) The combinatorial complexity of Vp 7 ^ (S) is Q(n 4 ) and 
0(n 4+E ) (for any e > 0). 

(b) If there is a unique closest pair p,q S S, then when 
c — > oo, the combinatorial complexity of Vn (S) is 
asymptotically 1. 

(c) For c > 0, the combinatorial complexity of Vp (S) is 
0{n 2+£ ) (for any e > 0). 

Proof: 

(a) To see the lower bound on the complexity of Vp 1 ] (5), 
note that every point on the segment pq has V-\- 
distance zero to the pair (p, q), and therefore, the inter- 
section of any pair of segments piqi and p2<?2 defined 
by sites Pi,qi,P2,Q2 G S is a feature of V p ] i (S). As 
is demonstrated in the proof of Theorem [8] the number 
of these features is £7(n 4 ). The upper bound is obtained 
by using the usual Davenport-Schinzel machinery, as in 
the proof of Theorem |2] 

(b) It is easy to verify that as c — > oo, the term c ■ \pq\ 
dominates the distance V c (v, (p,q)), and, hence, every 
point v in the plane is closer to the unique closest pair 
of sites p,q 6 S than to any other pair in S. Hence, the 
asymptotic diagram contains zero vertices, zero edges, 
and one face (the entire plane). 

(c) The proof is a generalized version of the proof of the 
special case c = 1. Refer to Figure [7] As in the proof 
of Theorem Q3] we assume that there is a point v in 
the region of (p,q), such that \vp\ < \vq\, and v is 
noncollinear with p and q. Our goal is to show that for 
any c > there exists a circle having q on its boundary 
and containing p, which is empty of any other site s, 
implying that p, q are Delaunay neighbors. 

As in the proof of Theorem Q3] let O v p\ be the locus 
of points q' for which V c (v,(p,q')) = V c (v,(p,q)). 
Thus, Ovpq is the Cartesian oval (v,p,c,k) consisting 
of all points q' that satisfy \vq'\ + c\pq'\ = k, where 
k = \vq\ + c\pq\ is constant. (Unless c = 1, this oval 
has exactly one axis of symmetry: the line joining the 
two foci v,p.) Then, if there were a site s within Oi c J„, 

fc) 

it would lead to a smaller value of V c , so O vpq must 
be empty of sites other than p. 

As before, let r be the ray emanating from q perpen- 

(c) 

dicular to and pointing into O vpq , and let o be the point 
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Figure 7. The Cartesian oval Ovpq is the locus of points q', for which 
\vq'\ + c • \pq'\ = \vq\ + c ■ \pq\. The ray r passes through q and is 
peipendiculai' to Ovpq , and intersects the axis of symmetry of O vpq at the 
point o. The circle C is centered at o, and is tangent to O v c J q at q. For any 
point x on the axis of abscissas residing inside C, t(x) denotes the point 
of C lying above x. 

where r crosses the line pv. 

Let us further suppose that c ^ 1. Without loss of 
generality, assume that 0% vq is symmetric with respect 
to the axis of abscissas (see Figure |7); consequently, 
the points p, v, and o belong to the latter. Let x p , x v , 
x , and x q denote the corresponding coordinate of p, 
v, o, and q, respectively. 

Consider a circle C centered at o of the radius R = \oq\. 
By construction, C is tangent to oi] q at q. 
For any x £ R, such that the point (x, 0) lies inside C, 
let t(x) denote the point of C lying above (x, 0). For 
any such x, let 

f v (x) = d(v,t(x)) 

= ^JB? - (x - x f + (x- x v y 
= y/2(x - x v ) ■ x + xl + x 2 v + R 2 . 

Since f v (x) represents a square root of a linear func- 
tion, it is concave on its domain. The same will hold 
for a function f p (x) = d(p,t(x)). Consequently, their 
weighted combination f(x) = f v (x) + c ■ f p (x) is also 
concave on the same domain, and, thus, has a single 
local maximum. 

Recall that the circle C is tangent to 0^ pq at q by 
construction. It is easy to see that C is tangent to 
Oy Pq from the inside: otherwise, x q would be a local 
minimum of f(x) achieved at an inner point of the 
domain, contradicting the concavity of f(x). It follows 
that f{x) has a local maximum at x q . Together with 
the previous observation, this implies that f(x) has a 
global maximum at x q . This means that q is the only 

fc) 

common point of Ov Pq and the upper half of C. By 
symmetry, we conclude that C lies inside oi c p \ and 
touches it at q and the point symmetric to q. Thus, C 
must be empty of sites other than p. 
It remains to demonstrate that p lies inside C. To 



this end, it is sufficient to show that the point o lies 
between v and p; then, as in the case of c = 1, the 
needed property can be easily derived using the triangle 
inequality. 

Let us argue as follows. The above reasoning can 
be carried out for any point q' G O c vpq noncollinear 
with v and p, providing us with a maximum empty 
circle inscribed in O c vpq , and tangent to it at precisely 
two points — namely, at q' and its symmetric point. 
It follows that the medial axis of 0% vq is a segment 
of the line Up through v and p. Let v' and p' be 

— (c) 

the intersection points of vp and Ov Pq being closer 
to v and p, respectively (see Figure |7). Consider the 
circle C v with radius \vv'\ centered at v. Obviously, 
v' is a common point of C v and oi c J q , but any other 
point z of C v lies strictly inside O c vpq , since for any 
such point z, we have \zv\ = \v'v\ and \zp\ < \v'p\. 
This implies that the radius of curvature of Ovp q at v' 
is greater than A similar statement holds for p' . 

Consequently, the two endpoints of the medial axis 
must lie between v and p, and the same must hold 
for the point o. 

We conclude that C is a circle containing both p 
and q and otherwise empty of sites, so p and q are 
Delaunay neighbors. Hence, there are Q(n) pairs of 
sites that generate regions in the Voronoi diagram, and 
the claim follows from the standard Davenport-Schinzel 
machinery. 

■ 

VIII. Conclusion 

In this paper, we have investigated 2-site Voronoi dia- 
grams of point sets with respect to a few geometric distance 
functions. The Voronoi structures obtained in this way 
cannot be explained in terms of the previously known kinds 
of Voronoi diagrams (which is the case for the 2-site distance 
functions thoroughly analyzed in J5]), what makes them 
particularly interesting. On the other hand, our results can 
be exploited to advance research on Voronoi diagram for 
segments. Potential directions for future work include con- 
sideration of other distance functions, and generalizations to 
higher dimensions and to fc-site Voronoi diagrams. 
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